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Abstract 

Repulsive gravity has its origin in the 1939 article of Oppenheimer and Sny- 
der which describes a collapsar, that is an idealized star of noninteracting mate- 
rial (dust) collapsing under its own gravity. The stellar material has a final state 
resembling a football, that is a significant part of it is concentrated in a thin 
surface shell. An interior pressure is exerted by the strong gravitational field, 
equivalent to a negative mass. However, the OS solution has been misunder- 
stood, the shell's position being incorrectly identified with the "event horizon" 
in black- hole theory. By expressing the OS metric in the preferred frame of 
harmonic coordinates, as advocated by N. Rosen and V. A. Fock, we obtain a 
description of the gravitational force field and thereby explain the football-like 
mass distribution. While half the material is concentrated in a shell occupying 
a small fraction of the radius, some material is spread throughout the interior, 
unlike the concentration in a black hole's singularity. We deal with the singu- 
larity in density at the shell surface by including Fermi pressure of degenerate 
electrons, for a shell density comparable to a solar mass-sized white dwarf. We 
show that use of the harmonic (inertial) metric but not the Schwarzschild metric 
leads to an OS-type solution. 



1 Introduction 



Repulsive gravity enables us to model supermassive galactic centres such as 
Sagittarius A* with most of the material concentrated in a shell which is inflated 
by a strong interior gravitational field. The mass density in the shell is similar 
to that in white dwarf stars, so we anticipate that degenerate electron Fermi 
pressure balances the repulsive gravitational field and spreads out the singularity 
in surface density of the dust star model. We call this a supermassive white- 
dwarf (SWD). Detailed investigation of SgrA* at the gravitational radius of 
OS, as is envisaged in the Event Horizon Telescope project, should discriminate 
very clearly between the SWD and black-hole models, as electromagnetic rays 
behave very differently in the two. Accretion of matter by SgrA* from the G2 
approaching dust cloud would also be different in the two models. 

The earliest indication that gravity may be repulsive is in the 1939 article of 
Oppenheimer and Snyder 1 (OS), describing an idealized collapsar composed of 
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" dust" , that is matter which distorts the gravitational metric but has no inter- 
action with itself. Penrose and subsequent experts in General Relativity (GR) 
failed to appreciate the novelty of the OS solution; it is in fact an immediate 
consequence of the OS metric that their collapsar has a stable end state whose 
radius is what they termed the gravitational radius - nowadays frequently re- 
ferred to as the Schwarzschild radius. It is only fair to point out that at least the 
better known of the OS authors should be given some part of the blame for the 
misunderstanding, since J. R. Oppenheimer had coauthored, with G. M. Volkoff 
a few months earlier an article^ (OV) which tried to analyse the collapse using 
Newtonian gravity forces. OS claimed that their results confirmed those of OV, 
whereas we show below that their consistently GR model in fact refutes it. 



2 The OS metric 



A rock, dug up 74 years ago, is wrongly classified as guartz and left on a museum 
shelf. After washing, it is reclassified as a huge diamond, about a hundred times 
the Koh-i-Noor. 

The surface of the OS collapsar describes a free-fall geodesic of the Schwarzschild 
metric, that is its radius ro (i) satisfies 

, 2 3/2 + 1 

t = --r -2^ + 1*1-— , (1) 



where we have chosen units in which the gravitational radius 2m is 1. The fact 
that t tends to plus infinity as ro tends to 1 led OS to infer at the end of their 
Abstract that "....an external observer sees the star asymptotically shrinking to 
its gravitational radius" . The interior r (t) < ro (t) is described by the metric 

ds^^<^- f^-^dR 2 - r 2 (d9 2 + sin 2 9d^) , (2) 

where R, with < R < 1, is a comoving coordinate, meaning that the interior 
free-fall geodesies are given by R ^constant; in particular R = 1 is identical 
with r = ro (t) . This latter condition, together with the requirement that the 
metric at R = 1 be continuous with the exterior Schwarzschild metric, led OS 
to put 

t{r,R) = - 2 -y^-2^ + \n^^ i , (3) 
where the cotime y is 

r R 2 - 1 

v{r > R)= R + ^2— • (4) 
Note that, while they referred to t as the "external time", this latter equation 
is an analytic continuation of t which establishes (r, t) as both exterior and 
interior coordinates. It is now a simple conclusion, which for some reason OS 
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failed to draw, that the stable end state, corresponding to t = +00, is obtained 
by putting y = 1, that is 

R(3-R 2 ) 

r °o = ^^ " ■ (5) 

Such a stable end state is in clear contradiction with the conclusion of Op- 
penheimer and Volkoff[3] and therefore with the opening sentence of the OS 
article itself. In contrast with black hole theory it gives a material density 
spread over the entire interior of the " event horizon" . They showed elsewhere 
in the same article that the initial state (t — > —00) has a uniform distribution 
of stellar material; this means, for example, that half of the total material is 
contained between R = 0.7937 and R = 1. Note that equation {3J informs us 
that r ~ R as t tends to minus infinity. Thus the OS collapsar has an end 
state with half of its material concentrated in a shell between r m = 0.9406 and 
foo — 1; we propose for it the title of the football collapsar. 

Actually the concentration of matter in the shell of the football is greater 
than this calculation indicates. The derivative of with respect to R contains 
a factor (1 — R), which indicates infinite density as R tends to 1. However, this is 
not as extreme as the singularity thought, by black-hole theorists, to lie at R = 0; 
rather, the physical conditions in the shell resemble those in a white dwarf of 
solar mass; a proper treatment of this region will require us to rehabilitate 
the notion of a gravitational force, and to incorporate both that force, which 
compresses stellar material into the shell, and the degenerate electron pressure 
of that material, which resists compression. That is a study which goes beyond 
the OS dust model, and which will require a different treatment. 

3 Modifying OS to meet a causality condition 



A small piece of quartz is separated from the diamond. 

We have just argued that the OS article does not support black-hole theory. 
However, the OS solution is problematic in that, as t approaches plus infinity, 
the collapsar reaches a certain stage when light cannot escape from it. OS gave 
minimal details to justify this anomaly, and we pointed out[4j that this feature 
of their model conflicts with the requirement of Hilbcrt causality enunciated in 
1917 by Hilbert[5], and revived more recently by Logunov and Mestvirishvili[6 . 
It requires only a simple correction to the OS metric [3], which leads to replacing 
©by 

r R 2 - 6R + 5 

V= R 4 ' (6) 

to ensure that radial light rays leaving R = for all t escape from the collapsar; 
this removes the problematic feature of the OS solution, which itself formed the 
basis for Penrose's black-hole theorem[2]. The modification also increases the 
shell concentration of the football, because 

-» = ^£ , (7) 
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which means that half of the material is now in 0.9659 < R < 1. 

Setting aside any predilection one might have for the geometrization of space- 
time, it is only natural to enquire "What pressure inflates the football?" In our 
dust collapsar there are no forces other than gravity, so once we have asked 
the forbidden question the only possible answer to it is that gravity itself, now 
turned from attraction to repulsion, inflates the football. A clear way to see 
repulsive gravity at work is to explore the trajectory of a test particle which 
crashes into the surface R — 1 at surface radius ro, with a speed greater than 
the speed at which the surface itself is contracting. We limit attention to radial 
trajectories for which the OS metric may be simplified, putting x — r/R, to 

ds 2 = xdx 2 - x 2 dR 2 , (8) 

for which there is an integral 

,dR 



x 
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= -C (C > 0) . (9) 



ds 

It then follows, putting xq = ro, that 

RM ^-L TO 4 ' • (10) 

and hence, putting y — 1 in ([6]), 

1 r x ° C 

Too = Xoo (3 - 2y / a^) , y/x^=l + - / —====dx . (11) 

1 J Xxi vx 6 + C z x 

We may thus compute the end point of the crash particle as a function of its 
time of arrival at the surface (effectively xq) and its crash speed (effectively C ). 
For the ultrarelativistic case, C — > +oo, the integration comes out in elementary 
functions and TVo takes the value 



= 1 (1 + V^) 2 (2 - v — i . 112) 



and, for all finite C, > r^, so that gives the maximum penetration of 
the football by a crash particle. Putting ro — xq — 1.45, which corresponds 
to the collapsar being only just inside its own photonsphere at r = 3m, gives 
r oo — 0.9666, and so, for all radii less than this, the crash particle fails to 
get beyond the "half-mass shell" which we obtained at the end of the previous 
paragraph. 

Thus, even ignoring collisions with the material of the football, such a crash 
particle encounters a sharp decelerating force on entering the shell, and that 
same force is what compresses the stellar material to form the shell. To make 
the gravitational force more concrete we shall now select a more appropriate 
coordinate frame. 
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4 Harmonic coordinates 



The facets of the diamond are cut and polished to reveal an object of wonderful 
simplicity. 

In general the cartesian harmonic coordinates ([7] pl63) x 1 ^ corresponding 
to the metric 

ds 2 = g^dx^dx" (13) 

are defined by 

□< = (^ = 0, 1,2, 3) , (14) 
where □ is the dalembertian operator 

n = -^=d^^— g d v . (15) 



-f) 



The Schwarzschild metric is 



ds 2 = rj -^-dt 2 s ~ -^-jdr 2 ~ r| (d9 2 + sin 2 6d<t> 2 ) , (16) 



and has the dalembertian 



r s -l Q2 2r s -l a , L 
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n = —?—dt g -- — # s -— V- d rs + — , (i7) 

rs - 1 s r s r s r s 

where 

L 2 = -d 2 g - cot <9<9 e - esc 2 Qd\ . (18) 
The spherical harmonic coordinates (i, r, 9, 4>), defined as those solutions of 

□t = 0, □ (rcos0)=O , (19) 

which arc asymptotically the same as the Schwarzschild coordinates for r$ — > 
+oo, are 

t = t s , r = r s -\ , (20) 

a result first obtained by Fock[5] and by Rosen[5] . Having designated the 
Schwarzschild metric with a subscript S, we may drop the subscript H of the 
harmonic metric, which is 

l9 2r-l , 2r + l , (2r+l) 2 , , 2 . . , 

ds 2 = -dt 2 - 7T—-dr 2 - y - — -i- (de 2 + sin 2 Odcj) 2 ) , (21) 

with the dalembertian 



2r-l d * 2r + l^ (2r + 1) 2 (2r + 1) 2 ' 
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It may be used to describe the gravitational field in the region r > r (t) , due 
to a spherical mass of radius r (t) ; the integral of the mass density over the 
interior of this sphere is m, which in our units takes the value 1/2 and is also 
the gravitational radius of the collapsar. Note that the latter quantity took the 
value 2m = 1 in the Schwarzschild metric so that, in the harmonic description 
as opposed to Schwarzschild's, the density of the collapsar is increased by a 
factor of 8. 

Now we pass to the OS metric. The substitution 

3 i? N 2 



converts it into 



r s = R[u+---) («>0) (23) 



ds 2 = ( u + - - [Adu 2 - 4dRdu - R 2 (d9 2 + sin 2 ddcj?)} (24) 



2 2 
with the dalembertian 
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□ 



(3 + 2u - R) 4 



(25) 

Just as R labels the particle geodesies, u labels the null geodesies, or outgoing 
light rays; this property of u is demonstrated by the absence of dR 2 va (EM)) , and 
of d 2 in □. The harmonic radius satisfies 

□ (rcos(9) = , (26) 

and, since it must be continuous with the exterior metric, r = r$ — 1/2 at R = 1. 
The solution is 

r = r s + - — = -- -(1- Rf + u(3R- R 2 ) +u 2 R (u > 0) .(27) 

Then the end state (u = 0) is 

roo = \-\{l-Rf , (28) 

for which half of the stellar material is in 0.4956 < < 0.5; in harmonic 
coordinates the density of the football shell is substantially greater than in 
Schwarzschild coordinates. The harmonic cotime y (R, v) must satisfy 

nt = o, i = -^3 /2 _ 2V y + ln ^±l (29) 

This gives the equation for y as 

y-3 fdy , dy\ dy =Q . (3Q) 



2y(y-l) \dR du J dR 
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where, at R = 1, both y itself and the metric (f2"4"|) must be continuous with the 
exterior metric, that is 



,2 



J( V A J2. „,2 ^2 , „ ; „2 flj^ 

y - 1 y 

which means that 



ds 2 = -^-dy 2 - -^-dr 2 - y 2 (d8 2 + sin 2 Odcj, 2 ) , (31) 
y - 1 y - 1 



y(l,y) = (y + l) 2 , J|(l |U ) = -u . (32) 

It should be noted that the causality condition we imposed in [3], namely 
that y (R, 0) = 1 for < R < 1, ensures that the solution of (|30l) is nonsingular 
at u = 0. In the original OS model a similar conversion to harmonic coordinates 
is not possible; Hilbert causality is essential for such a conversion. Equation (|30[) 
may be integrated numerically, starting from the asymptotic solution as u tends 
to 0, namely 

r i + JL (2ar w, a{X )= * 2 ( 3 + x ) (w ^o,o<x<i) , 

(1-X)(2 + X) 2 v ' 

(33) 

where we have put R = 1 — X. The corresponding expression of r, for all u, is 
r = 1 ~ X + u (2 - X - X 2 ) +u 2 (l-X) (u>0,0<X < 1) . (34) 



5 A supermassive white dwarf 

White dwarfs of solar mass have been understood from the 1930s in terms of 
the Chandrasekhar|:13 analysis, which (see, for example, [7] (11.3)) balances 
the Fermi pressure of a degenerate electron gas against the gravitational force. 
Although gravity in a white dwarf is very strong, and may indeed eclipse the 
Fermi pressure, it nevertheless remains Newtonian. The Oppenheimer articles 
of 1939 [3 [1] were two early attempts to extend that analysis to more radical col- 
lapse situations, but now, many decades later, nobody has yet managed to treat 
material forces, like the Fermi pressure, in the context of Einstein's gravitation 
theory. 

A certain generalization may be achieved ([7] (11.1)) by considering the GR 
field equations with the metric 

ds 2 = B(r)dt 2 - A(r)dr 2 - C(r){d9 2 + sin 2 0d(/) 2 ) , (35) 
and the stress tensor for a perfect fluid 

T ^ = pgs* + ( p + p ) u» U v , (36) 
with its associated equation of hydrostatic equilibrium when U = (1, 0, 0, 0) 

B' -2p' , s 

= — - ■ 37 

B p + p 
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This latter equation, together with the two field equations 
and 

Rss + lA R - + ^ R - = - 1 + Tr^A=- 87rGpC > (39) 
and a suitable equation of state p = p (p, r) , provides us with four equations for 
five unknowns, so we are lacking one equation. 

There is a well known procedure for simplifying the metric, using the freedom 
of coordinate choice in GR, so that C (r) = r 2 . In that case A (r) and B (r) may 
be eliminated from the above three field equations to give ([7] (11.1.12) and 
(11.1.13)) 

d r 2 p' + AirGr^p (p + p) 



dr p + p — 2rp' 
In the Newtonian regime this becomes 



AnGr 2 p . (40) 



d r2? ' = -^Gr 2 p , (41) 



dr p 

and an analytic solution of this latter equation is possible for a zero-temperature 
dilute^ electron gas, having a polytropic equation of state with index 5/3. The 
solution in that case is obtained by integrating from r = 0, with the initial 
conditions 

p(0)=p , p'(0)=0 , (42) 

and shows p monotone decreasing to a zero at n = r% (p ) . If we were to attempt 
to use the same equation of state in (|40|) , we would arrive at a similar conclusion; 
the density decreases monotonically from p at r = to zero at n = ri (Po)- So 
there is nothing resembling the final state of the football collapsar. However, 
it would be an error to combine (|40| with the polytropic equation of state, 
because the latter refers to a specific inertial frame of coordinates, and there is 
no reason to presume that the interior Schwarzschild metric defined by C = r 2 
is that frame. 

Now we propose to investigate what happens if instead we use the harmonic 
coordinate r in conjunction with the polytrope equation of state, so that C (r) 
is given in terms of A(r) and B{r) by ([7] (8.1.15) ) 

i{cM\=*nTAB . (43) 



dr V V A 
We follow Weinberg ([7] (11.3)) and define 



1/2 h 2 /^2 n \ 2/3 

(44) 



onGpi J Z()m e mr)iT \ mo 



1,1 Dilute" means the volume per electron is large compared with the cube of the Compton 
wave length h/mc. 
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where m e is the electron mass and mo is the mass of nuclear material per elec- 
tron, which is approximately one deuteron. However, the calibrating quantity 
p x is now defined as the maximum value of the density rather than the density 
at the centre. Then the polytrope equation becomes 



p = F9p, F 



Pi 



5m D \p c 



2/3 



(45) 



where p c , which is a critical density at which the energy of an electron at the 
Fermi surface equals y/2 times its rest mass, that is 



Pc 



&7TUI 



d ^m e cy 



The equation of equilibrium now becomes 

5F d9 



l_dB 
B~d7 



l + F9dr 



(46) 



(47) 



Note that the parameter F does not reach order 1 until p 1 exceeds p c by a 
factor of order 10 6 , that is long after its polytrope equation has changed index 
from 5/3 to 4/3. 
Now, defining 



(48) 



and using a prime henceforward to denote differentiation with respect to £, the 
harmonic condition becomes 



d_ 



J-(l + F6)- 5 



and the two field equations become 



^ (1 + Fey 5 = 5F 7 (1 + F9) 9 3/2 



and 



d_ 



1 - 5F 7 6» 3/2 



(49) 

(50) 
(51) 



These three equations may be written as 

(l + ^6>) 6 6' 3/2 





10Fa 2 7 


a' = 




7' 




1 + F9 


9' = 




5F 


7" = 


2a (1 + j 



a' 
2a~ 



a 
2a~ 



(52) 
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where 

13 = y - 2a£ (1 + F6f , (53) 

and in this form may be integrated numerically with the initial conditions 9 (0) = 
7 (0) = 7' (0) = 0, and a (0) treated for the moment, like F, as an adjustable 
parameter. Exactly as in the Newtonian theory, we stop when 9 (£) reaches its 
first zero ^ , which, as in that theory, is a function of p l through the parameter 
F, but also now of ao- To complete the solution ao must now be adjusted so 
that the metric is continuous with the exterior harmonic metric at £ = £1, that 
is 

A(r 1 )=B(r 1 )- l = ^^ ) C(r 1 ) = (r 1 +mf , (54) 
T\ — m 

where m is the integrated mass from r = to r\. In terms of the rescaled 
variables a (£) and 7 (£) , since B is determined up to a constant factor by 9 (£), 
these conditions reduce to the single one 



afo)= (55) 

and may be met by adjusting a(0), which must be greater than 1 in order to 
ensure that 7 > £ 2 in the entire interval < £ < The solution then depends 
on the single parameter F, that is p 1 , to which m and r\ are related as in the 
Newtonian theory. 

Our conclusion is that assuming a polytrope with index 5/3 for the equa- 
tion of state in harmonic coordinates must lead to a metric showing repulsive 
gravity at the centre; this follows directly from ([4"?| . because now 9' (0) is pos- 
itive which means B' (0) is negative. On the other hand we saw above that in 
Schwarzschild coordinates the polytrope equation of state gave a density mono- 
tonically decreasing from the centre to the surface. 



6 Conclusion 

We saw in sections 2 and 3 that in a dustball collapsar the density decreases 
towards its centre. Section 5 establishes that same behaviour when an elec- 
tron gas is incorporated in the model. This validates our supermassive white 
dwarf (SWD) concept. The SWD has relativistic gravity as distinct from a nor- 
mal white dwarf with Newtonian gravity ([17], Investigation 12.4) and density 
~1 tonne/cm"3. Note that the higher density white-dwarf regime ([17], Inves- 
tigation 12.5) has polytropic index 4/3 rather than the 5/3 used here. The 
Oppenheimer-Volkoff modelling [3] also assumed Newtonian gravity, so cannot 
apply to our SWD or to black hole analysis. 

We found the OS analysis had to be modified to satisfy (Hilbert) causality, 
ensuring connection by light rays is maintained throughout the collapse. We 
also found that the Schwarzschild coordinate frame gave an unphysical result; 
since it is one of an infinite set of frames permitted in GR, we can with some 
confidence assert that the polytropic equation of state is almost certainly not 
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valid in it. We instead selected the harmonic frame of section 4, imposing 
the harmonic condition on C (r) given by Weinberg (8.1.15). We have stated 
the case for this choice in previous articles QZ3 [H] [H] ; our motivations may be 
grouped under its computational simplicity [7], its historic basis, especially in its 
use by Einstein to found the study of gravitational waves [13). and its status as 
a possible cornerstone for a field theory of gravitation[9"] |15) |16j . 

The radio and infrared emissions from SgrA* are thought to emanate from 
gas and dust heated to millions of degrees while falling under the central gravity. 
In the SWD model such heating would occur just above the shell surface in a 
corona. The corona would be created by the inputs of dust, including gas 
evaporated from it, and depend on emission-cooling; these complexities could be 
modelled, whereas the alternative "interaction with a massive source of gravity" 
being source of the observed emissions is vague and speculative [18] 

The X-ray flaring of SgrA* and other distant galaxy centres has been hy- 
pothesised as arising from the shredding and consuming of a passing star [115] • 
In the SWD model, this would need a large concentration of matter to impact 
the shell and temporarily destabilise it - a stellar mass of WD density - as we 
estimate a solar-sized impactor would be a factor 10 000 times less massive than 
the impacted shell. On the other hand, the expected impact in 2013 on SgrA* of 
the G2 dust/gas cloud of a few Earth masses would be far too small to perturb 
the football shell. With a closest approach some 3000 times the gravitational 
radius (0.1 AU) and dimension ~1AU, it could hardly affect even a large corona 
of size "0.1AU. If SgrA* possesses a much larger accretion disc, this could ab- 
sorb material from the G2 cloud. However, the accretion onto the SWD from 
such an accretion disc would be over a longer time-scale and would not give the 
predicted fireworks display. 
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